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Abstract 

We present a summary update of the QCD spectral sum rule (QSSR) results for the running and 
perturbative pole quark masses, the fo and fs leptonic decay constants, the heavy-to-light and 
heavy-to-heavy exclusive transition-form factors. Analytic expressions of these latter quantities 
are presented, which give a deeper understanding of their g^- and infinite mass-behaviours. A 
short comparison of the QSSR results with alternative approaches is done. 
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1 QCD Spectral Sum Rules (QSSR) 



QCD spectral sum rule (QSSR) a la SVZ ||T| (for a recent review, see e.g. 0) has shown since 
15 years, its impressive ability for describing the complex phenomena of hadronic physics with 
the few universal "fundamental" parameters of the QCD Lagrangian (QCD coupling a^, quark 
masses and vacuum condensates built from the quarks and/or gluon fields), without waiting 
for a complete understanding of the confinement problem. In the example of the two-point 
correlator: 

n,(g2) ^ I J d\ e^^^ (0|rj,(x) (J,(o))^ |0), (1) 

associated to the generic hadronic current: Jfe(x) = qVh{x) of the q and 6-quarks (F is a Dirac 
matrix which specifies the hadron quantum numbers), the SVZ-expansion reads: 

D=o,2,...dimo=D [-lyh ~Q ) 

where /i is an arbitrary scale that separates the long- and short-distance dynamics; C^'^-' are the 
Wilson coefficients calculable in perturbative QCD by means of Feynman diagrams techniques; 
(O) are the non-perturbative condensates of dimension D built from the quarks or/and gluon 
fields [D = corresponds to the case of the naive perturbative contribution). Owing to gauge 
invariance, the lowest dimension condensates that can be formed are the D = A light quark 
mqi^pip) and gluon {agC^) ones, where the former is fixed by the pion PCAC relation, whilst the 
latter is known to be (0.07 ± 0.01) GeV^ from more recent analysis of the light and heavy 
quark systems 0. The validity of the SVZ-postulate has been understood formally, using 
renormalon techniques (absorption of the IR renormalon ambiguity into the definitions of the 
condensates, UV renormalon cannot induce some extra 1/M^-terms not included in the OPE) 
0, and/or by building renormalization-invariant combinations of the condensates (Appendix 
of 10 and references therein). The SVZ expansion is phenomenologically confirmed from the 
unexpected accurate determination of the QCD coupling a.^ and from a measurement of the 
condensates from semi-inclusive tau decays 0. The previous QCD information is transmitted 
to the data through the spectral function lmlliy{t) via the Kallen-Lehmann dispersion relation 
{global duality) obeyed by the hadronic correlators, which can be improved from the uses of 
either a finite number of derivatives and finite values of (moment sum rules): 



f°° dt I ^ ^ , ^ 

/ -imn,(t), 3 



or an infinite number of derivatives and infinite values of g^, but keeping their ratio fixed as 
r = n/q^ (Laplace or exponential sum rules): 

rco 1 

C(t,M^) = / dt exp(-tr) - Imnf,(t), (4) 

JM^ TT 

for rriq = 0. Non-relativistic versions of these two sum rules are convenient quantities to work 
with, in the large-quark-mass limit, after introducing the non-relativistic variables E and tn'- 
t = (E + Mbf and ttv ~ MbT. In the previous sum rules, the weight factors enhance the contri- 
bution of the lowest ground-state meson to the spectral integral, such that, the simple duality 
ansatz parametrization: "one narrow resonance" + ^^QCD continuum" , from a threshold tc, 
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gives a very good description of the spectral integral. The previous naive parametrization has 
been tested successfully in the light-quark channel from the e+e~ — > / = 1 hadron data and 
in the heavy-quark ones from the e^e~ ^ ip or T data, within a good accuracy. In princi- 
ple, the pairs {n,tc), {T,tc) are free external parameters in the analysis, so that the optimal 
result should be insensitive to their variations. Stability criteria, which are equivalent to the 
variational method, state that the best results should be obtained at the minimas or at the 
inflexion points in n or r, while stability in tc is useful to control the sensitivity of the result 
in the changes of tc-values. To these stability criteria can be added constraints from local du- 
ality FESR, which correlate the tc- values to those of the ground state mass and coupling p. 
Stability criteria have also been tested in models such as the harmonic oscillator [^j , where the 
exact and approximate solutions are known. However, though / would personally expect that 
the true result is obtained near the beginning of the tc-stability region [| despite the fact that 
in some cases this value of tc is larger than the phenomenological guessed position of the next 
radial excitation, one can fairly state that the most conservative optimization criteria, which 
include various types of optimizations in the literature, are the one obtained in the region, 
starting from the beginning of r/n stability 0, until the beginning of the tc stability. One can 
a posteriori check that, at the stability point, where we have an equilibrium between the con- 
tinuum and non-perturbative contributions, which are both small, the OPE is still convergent 
and the expansion certainly makes sense. The results which will be quoted below have been 
obtained within the previous stability criteria. 



2 The heavy- quark- mass values 

Many efforts have been devoted to the study of the quark masses Using the present 



world average value as{Mz) = 0.118 ±0.006 ITTI, the first direct determination of the nmning 



mass to two loops, from the \l/ and T systems, is [|T2| 



Me(Mf^2) = (1.23+°:°^ ± 0.03) GeV 

m,{Mr') = (4.231^4 ± 0.02) GeV, (5) 
where the errors are respectively due to and to the gluon condensate. Using the relation 

mi: 

M,=^QiMi)[l + l{^)+K,[^yy (6) 

where Kh ~ 12.4, Kc — 13.3 [|1^, one can transform this result into the perturbative pole mass 
and obtain, to two-: 

Mf = (1.42 ± 0.03) GeV 

^2 _ (4_g2 ± 0.02) GeV , (7) 



^ This value of tc is about the one fixed by FESR duahty constraints. In this region, the result is certainly 
insensitive to the form of the continuum model. 

^This corresponds in most of the cases to the so-called plateau often discussed in the literature, but in my 
opinion, the interpretation of this plateau as a sign of a good continuum model is not sufficient, in the sense 
that the flatness of the curve extends in the uninteresting high-energy region where the properties of the ground 
state are lost. 
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and three-loop accuracy: 



Mi 



(1.62 ±0.07 ±0.03) GeV 
(4.87 ±0.05 ±0.02) GeV. 



It is informative to compare these values with the ones of the pole masses from HQET |jT5|, ^ 



and non-relativistic sum rules to two loops [O: 



M, 
M, 



NR 



NR 



:i.45l°:°^±0.03) GeV 
± 0.02) GeV, 



(4.69 



(9) 



where one might interpret the small mass difference of about 70 MeV as the size of the renor- 
malon effect into the pole mass. Indeed, an explicit resummation of the leading (/Ja^)" terms 
increases the previous two-loop estimate by about 100~200 MeV It also indicates that the 
three-loop value in Eq. (8) already gives a good estimate of the pole mass to all orders of PT. 
Eq. (8) compares quite well with the dressed mass M^'' = (4.94 ± 0.10 ± 0.03) GeV, obtained 
from a non-relativistic higher-order Balmer formula based on a bb Coulomb potential [|17| and 
with lattice calculations |l^. One can also use the previous results, in order to deduce the 
(non)-relativistic pole mass-difference of the b- and c-quarks both evaluated at M^: 



Mb-M, 



--Mi 



(3.54 ± 0.05) GeV. 



(10) 



Finally, the lesson which we can learn from the previous discussion is that one should be very 
careful in using the numerical value of the quark mass. Indeed, for consistency, one should 
first understand the definition of the mass used in the analysis and know to what loop-accuracy 
the analysis is done (for a recent compilation of the running quark masses, see e.g. |T^.). 



3 The decay constants and the B^-parameter 

The decay constants fp of a pseudoscalar meson P are defined as: 

(m, ± MQ)(0|g(275)Q|P) = V2Mlfp, (11) 

where in this normalization /^^ = 93.3 MeV. A lot of efforts have also been done for the estimate 
of these decay constants [0]. However, the most pertinent result has been obtained, for the first 
time, in pO[| : 

Jd^Jb^ 1.4/. , (12) 

which clearly shows a violation of the Infinite Mass Effective Theory (IMET) 1/ a/ Mfe-scaling 
law. Several obscure and unjustified criticisms have been adressed later on in order to discredit 
such a result, but a numerical estimate of the 1/Mfe-corrections complemented by the lattice 
results and by the analytic estimate from HQET and semi-local duality has provided a better 
understanding of the unexpected result in Eq. (12). Indeed, using the present update best 



estimate from the Laplace sum rule pT| : 



fo ^ (1.35 ±0.04 ±0.06)/. 

/b ^ (1.49 ± 0.06 ± 0.05)/. , (13) 
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consistent with the value of the relativistic two-loop pole mass given previously, and using the 
value of the decay constant in the static limit |16|, |22| : 

/^^ (1.98 ±0.31)/^, (14) 

the l/Mfe-corrections are found to be [2^, |2^, |l^, p2| : 



fs^^b^ (0.33±0.06)GeV^/2ay^i|l - 
2 a, (A -1.1 GeV) {B - 0.7 GeV^) | 

3V ^ M| / ^^^^ 

which, one can qualitatively compare with the one obtained from the analytic expression of 



the moments or from the semilocal duality sum rule leading to the interpolating formula : 



2a 3 El Tc'iuu) \ 



3 TT 88 M2 2 



and gives for Ec — 1.3 GeV: 



A ^ ^(Ms-M6)~lGeV, 

B ^ lEl-^iMB-Mb)'^0.5GeV', (17) 

The SU (3)-breaking effects to these decay constants have also been estimated analytically to 
be 0: 

^~^~(1.16±0.04)A, (18) 
Id Jb 

which is in good agreement with the range of values obtained from different lattice groups . 
The corresponding value of fn^ ~ (1.55±0.10)/,r is still compatible (within the errors) with the 
recent (indirect) measurements from WA75 and CLEOII (see e.g. |^), which need to be tested 
from direct (though difficult) measurements of the leptonic widths. However, independently 
of the charm quark mass-value, which affects strongly the value of fo Ull (/d increases for 



decreasing Mc), a value of fo^ larger than the rigorous upper bound of 2.14/^ deduced from 



28| and Eq. (18) is unlikely from the QSSR approach within the standard SVZ-expansion. 
We have also tested the validity of the vacuum saturation for the -B^-parameter, and we found 
that the radiative corrections for the non-factorized correlators are quite small (less than 15%), 



from which we deduce 29 



5b ~1 ±0.15. (19) 



4 Heavy-to-light transition-form factors 

One can extend the analysis done for the two-point correlator to the more complicated case 
of three-point function, in order to study the form factors related to the B Tx{p)lv and 
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B K*'~) rare decays. In so doing, one can consider the generic three-point function: 

V=-jd\ dS e*(^'^-^^) (0|TJi(x)O(0)4(i/)|0), (20) 

where Jl, Jb are the currents of the hght and B mesons; O is the weak operator specific 
for each process (penguin for the -ft'*7, weak current for the semileptonic) ; g = p — p' is the 
momentum transfer. The vertex obeys a double dispersion relation, which can be improved 
into the unique hybrid sum rule (HSR) |Q ^ : 

n{n, r') = - l-^ ds' ^ Im\/(., .'), (21) 



b 



corresponding to a finite number n (ra ^ 1 — 2 in the present processes) of the moments for 
the heavy-quark channel and to the Laplace for the light one. We have studied analytically 
the different form factors Q entering the previous processes [^, and we found that they are 



dominated universally, for Mi, oo, by the the light-quark-condensate contribution as: 

f(0)~^|l + 5i|. (22) 



where Xp is the integral from the perturbative triangle graph, which is constant as t'^Ed (dd) 
{t'^ and Ec are the continuum thresholds of the light and h quarks) for large values of M^. Unlike 
the case of /s, where the perturbative graph and the (gg) condensate are of the same order in 
Mb, the present dominance of the {qq) condensate allows a good separation of the lowest ground 
state contribution from the radial excitation. It also indicates that at = and to leading 
order in 1/M;,, all form factors behave like a/M^, although, in most cases, the coefficients of 
l/Mft due mainly to fs and of 1/M^ due to the perturbative graph are large Q. In the particular 
case of -B ^ vr/z/, the form factor can be simply written to leading order p3|: 




/+(0)^ 71^ 4^ U 0.15, (23) 



to be compared with the numerical estimate 0.25 (the factor reflects the off-shellness of the 
pion) and to the pion coupling to hadron pairs |@] . The study of the g^ behaviours of the form 
factors shows that, with the exception of the Ai form factor, their g^ dependence is only due 
to the non-leading perturbative graph, so that for Mb oo, these form factors remain 

constant from g^ = to g^^^. and have a weaker g^-dependence ( polynomial in g^ that can 
be resummed), than the pole model at finite M^ (here the value of pole mass which fits the 
form factors is about 5-6 GeV (see also |^)). The resulting Mj, behaviour at g^^^ is the 



one expected from the heavy quark symmetry. The situation for the Ai is drastically different 
from the other ones and from the pole parametrization. Here the Wilson coefficient of the (dd) 
condensate contains a g^ dependence with a wrong sign and reads: 



(24) 



■^The popular double exponential sum rule is not appropriate here as in this sum rule the OPE blows up for 
Mb — > oo. 

^In the standard notations, the relevant form factors are /+, {Ai, A2, V) for B ■n{p)lv and Fi for 
B —f K*j decays. 

^This feature also indicates that it is dangerous to extrapolate the Mh-dependence obtained at the c-quark 
mass to higher quark mass values. 
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which, for g^^^ = {Mb — Mp)'^, gives the expected l/^/Mb behaviour. This result also explains 
the numerical observation in ||3^ (similar conclusions using alternative approaches have also 
been reached in |^). Numerically, we obtain at = 0, the value of the B p{K*)'^ form 
factors: 

F]^""'' ^ 0.27 ± 0.03, ^^^^ -1.14 ±0.02, (25) 

which leads to the branching ratio (4.5 ± 1.1) x 10~^, in perfect agreement with the CLEO data, 
while the numerical agreement with the estimate in ||3^ from light cone sum rule (for a criticism 
on the unreliability for the construction of the hadronic wave functions on the light-cone, see 
e.g. 1^) may only be due to the importance of the perturbative contribution at this scale 
^ For the semileptonic decays, a determination of the ratios of the form factors gives a more 
precise prediction than from a direct estimate of the absolute values: 

44^. ^.1.11 ±0.01, 

v4i(0) Ai(0) 

^i^^ ~ 1.18 ±0.06, 444 - 1-40 ±0.06. (26) 
F^-^iS)) /+(0) ^ ' 

Combining these results with the "world average" value of /+(0) = 0.25 ± 0.02 and the one of 
Fi~''^{0), one can deduce the rate and polarization: 

~ (4.3 ± 0.7)|K6p X 10^2 g-i ^ 0.9 ± 0.2 

r r 

0.20 ±0.01 a = 2-^ - 1 ^ -0.6. (27) 

1 _ 1 7^ 

These precise results may indicate that, Vub can be reached with a good accuracy from the 
exclusive modes. The non-pole behaviour of Ai affects strongly the different estimates in Eq. 
(27), in particular the ones of Fp/F^ and a, such that a firm prediction of these quantities needs 
a an improved good control of the g^-dependence of the corresponding form factors. We extend 
the previous analysis to the estimate of the SU{3) breaking in the ratio of the form factors: 

Rp ^ /;^^(0)//;-(0), (28) 

where P = B, D. Its analytic expression is given in ||3^ and leads to the numerical result: 

= 1.007 ±0.020 i?D = 1-102 ±0.007, (29) 

which is typically of the same size as the one of fo^ and of the B K*-^ discussed before and 
reinforces the credibility of the present estimate. However, it is quite surprising that using the 
previous value of Rp, into the present value of the CLEO data ||40[| : 

one deduces: 

Vcd/V,s = 0.322 ± 0.056, (31) 

which is much larger than the value 0.226 ± 0.005 derived from the unitarity of the CKM 
matrix. This apparent discrepancy needs a further measurement of the previous process before 
a firm conclusion can be drawn (recall that MARKIII data |^1]] would imply a value 0.25 ±0.15 
compatible within the errors with Eq. (31).). 



^Lattice results on the radiative and semi-leptonic decays are reported in |39 
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5 B*Bn{'y) couplings and D* Dn{'y) widths 



As has been studied recently in , the previous processes are very similar to the other heavy- 
to-heavy transitions as they are dominated by the perturbative graph contributions. The non- 
leading 1/Mb corrections for the radiative decays are large as they come mainly from the heavy 
quark component of the electromagnetic current. This contribution is essential for explaining 
the large charge dependence in the observed radiative decay widths. For the B* meson, our 
predictions without any free parameters are: 

gB*B. ^ 14 ± 5 TB*-.B-j/rB'-^B-^ ^ 2.5, (32) 

where the latter indicates a large isospin violation, which deviates strongly from the naive 
static limit e^/e^ expectation, therefore showing the importance of the l/Mb corrections in this 
channel. For the D*-one, we find: 



while: 



,D0^- ~ 1.541^,0^^0^0 ~ (8 ± 5) keV, (33) 



Tn^-^n-y ^ (0.09+°:^°) keV 

r^.o^^o^ ^ (3.7± 1.2) keV. (34) 

The resulting total widths TD*-^aU — (12 ± 7) keV and TD*o^aU — (H ± 4) keV are much 
smaller than the present experimental upper limits. Improved measurements of these widths in 
the next r-charm factory machine should provide a decisive test of the predictions given here 
and should also help to clarify the disagreements among the present theoretical predictions. 



6 Slope of the Isgur— Wise function and Vcb 

From the QSSR expression of the universal Isgur-Wise function, to leading order in I /Mb [0: 

W(^-W) = (^y|i + ^^/(,) 

- {ddyg{y) + {a,G^yh{y)+g{dGd)T'k{y)'^, (35) 

where r is the Laplace sum rule variable and /, h and k are analytic functions of y. From this 
expression, one can derive the analytic form of the slope of the IW-function |H3| : 



Cphysiy = 1) ~ -1 + 6,ert + SnP ^ -I ± 0.02, (36) 

where at the r-stability region: 

Spert ^ -Snp ^ -0.04, (37) 

which shows the near-cancellation of the non-leading corrections, and we have added a generous 
50% error of 0.02 for the correction terms. This result is in agreement with the improved bound 



of Taron-de Rafael on the slope of the form factor |[44|| : 



F'{vv' = 1) > -1.5, (38) 
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based on the analyticity, the positivity and a mapping technology of the elastic 6-number form 
factor F defined as: 

{B{p')\h^b\B{b)) = {p + p'mq'). (39) 

and normalized as -F(O) = 1 in the large mass limit Mb — Md- The inclusion of the effects of 
the T states below BB thresholds by using the sum of the TBB couplings of 0.34 ± 0.02 from 
QSSR improves slightly this bound to: 

F'{vv' = 1) > -1.34. (40) 

Using the relation of the form factor with the slope of the Isgur-Wise function, which differs 
by —16/75 log as(M;,) one can deduce the final bound: 

C'(l) > -1.04. (41) 

The QSSR results are in good agreement with some other significant estimates given in the 
literature. 

Let us now discuss the effects due to the 1/M corrections, which can be done in two ways: By 
calculating the 1/Mf, corrections within HQET, (resp. [^) obtains: 

ri^((l) = 0.93 ± 0.03 (resp 0.89 ± 0.03) . (42) 

Here, the model-dependence enters when extrapolating the data at y = 1, and leads to: 

Vcb ^ (39.9 ± 2.9) X 10-3 (resp (41.7 ± 3) x lO'^), (43) 

Alternatively, one can use the value of the form factor at = (y = 1.5) from the sum rule in 
the full theory : 

F(1.5) = 0.53 ±0.09 (44) 
and the data in the whole range of y in order to deduce the slope: 

p2 = -C ^ 0.76 ± 0.2, (45) 

from a linear parametrization of the form factor. The model dependence enters in this analysis 
through the curvature of the form factor. The main advantage of this approach is that it does 
not rely on the previous theoretical conflict at y = 1. Using the different data from CLEOII, 
ARGUS and ALEPH @ , we obtain the average: 



'cb 



(39.9 ± 1.2 ± 1.4) X 10-^ (46) 



where the first error is from the data, and the second one is from the type of model-parametrizations. 
This result is in good agreement with the one in Eq. (43) and from the inclusive decays. 



7 Conclusion 



We have shortly presented different results from QCD spectral sum rules in the heavy-quark 
sector, which are useful for further theoretical studies of the 5-physics and which complement 
the results from alternative non-perturbative approaches. From the experimental point of view, 
QSSR predictions agree with available data, but they also lead to some new features, which 
need to be tested in forthcoming experiments. 
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